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We present two of our efforts directed toward the numerical analysis of neutron star merg- 
ers, which are the most plausible sources for gravitational wave detectors that should 
begin operating in the near future. First we present Newtonian 3D simulations including 
radiation reaction (2.5PN) effects. We discuss the gravitational wave signals and lumi- 
nosity from the merger with/without radiation reaction effects. Second we present the 
matching problem between post-Newtonian formulations and general relativity in numer- 
ical treatments. We prepare a spherical, static neutron star in a post-Newtonian matched 
spacetime, and find that discontinuities at the matching surface become smoothed out 
during fully relativistic evolution if we use a proper slicing condition. 



1 Introduction 

Mergers of neutron stars in binary systems are considered excellent astrophysical 
laboratories for gravitational wave astronomy, nuclear astrophysics and relativistic 
astrophysics. The gravitational waves emitted from these mergers are expected to 
be observed by gravitational wave detectors coming on-line in the next decade, such 
as LIGO, VIRGO, GEO and TAMA. Detailed predictions of waveforms are desired 
and need to be produced analytically and numerically for the extraction of physical 
information. 

Post-Newtonian (PN) calculations! 1 ] can accurately describe the evolution of the 
system when the separation between the two stars in the binary system is much 
larger than the stellar radii. The final phase of the neutron star mergers (NSM) 
requires fully general relativistic (GR) numerical simulations. In this paper, we 
present two of our efforts at building a bridge between these two approaches. 



2 Newtonian simulation of NSM including radiation reaction 

First we present,Newtonian 3-dimensional simulations including-.radiation reaction 
(2.5PN) effects cl. We note that the Kyoto and Max-Planck □ groups have also 
presented results of their NSM simulations using Newtonian equations plus 2.5PN 
corrections. We use the same set of equations, whicluare a reduction of the PN- 
hydro formulation of Blanchet, Damour and SchaferB. We are studying various 
initial conditions and comparing the effects of different equations of state. 

We evolve the Euler equations using a modification of the ZEUS-2D algorithm, 
staggered grid structure, second order van Leer monotonic interpolator, and Nor- 
man's consistent advection method a. We have conducted convergence studies to 
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delineate how spatial and temporal resolution affect the conservation of angular mo- 
mentum and energy in these models. In order to get sufficient angular momentum 
conservation, we find that we must update the "advection" terms before solving the 
Poisson equations. To solve the Poisson equations, we use a full multi-grid W-cycle 
algorithm. 

Here, we show an example of the gravitational waveform emitted by the coa- 
lescence of two equal mass stars. The equation of state is polytropic P = Kp 2 with 
initial mass 1.4M Q and radii i?* = 9.56Km. The stars are separeted initially by 
2.9i?*, are nonrotating, and have an orbital velocity taken to be the Kepler velocity 
with infalling radial velocity determined by radiation reaction This simulation was 
made with 129 3 grid zones (Ax — 0.893Km) and required 160 hrs of CPU time by 
Origin 2000 with an average time step 1.43 x 10 _3 ms for a total physical time of 
5.40 ms. 

In Fig.l, we show the gravitational waveform h x and luminosity with and with- 
out radiation reaction term. We see that the merger occurs earlier when we include 
radiation reaction terms. This may be explained by the loss of angular momentum 
frorn_the system induced by radiation reaction. Detail reports are now in prepara- 
tion!. 
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Fig. 1: Waveform h x and luminocity. Solid line is Newtonian and 
dashed line is with radiation reaction. 

3 GR evolution of a neutron star with PN matched initial data 

Our second effort in linking PN calculations to full GR numerical simulations in- 
volves construction of PN initial data for fully relativistic numerical treatment. 

As a first step, we constructed a single spherically symmetric neutron star by 
solving the Tolman-Oppenheimer-Volkoff equations; this solution is matched to a 
1PN vacuum metric at some surface outside the star. We then evolved this initial 
data using a fully relativistic spherical code, and investigated what happens at the 
matching surface. 

In the central region of the star, we set our fully relativistic metric to 

ds 2 = _ e 2*(R) dt 2 + e 2A(fl) dr 2 + R 2 {r) ^2 + ^2 ^2) _ (1) 
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where r = r(R). Under the harmonic gauge condition 4- [e* A p 2 ] — 2re* +A , we 
solve the hydrostatic equilibrium (TOV) equation and field equations assuming a 
polytropic equation of state. We define a "mass" M, given by M — Air f^ ' p(R)R 2 dR, 
where i?» is the radius of the star. We set the exterior metric of the neutron star 
by $(i?) = \ln{\ -^jf) + C, and e 2A = (1 - ^f-)~\ where R = r + M and the 
constant term C is to be fixed by matching later. The outer metric (PN region; 
r > a > where a is the matching radius) is given by 

ds 2 = -(1-2*7 + 2U 2 + [3U 3 )dt 2 + [1 + 2U + 2U 2 + \aU 3 ]dr 2 

+ [1 + 2U + U 2 - \aU 3 ]r 2 (d9 2 + sm6d<p 2 ) , (2) 

where U = rh/r, with rh the Kepler- measured mass of the star, [3 the coefficient 
of a 2PN term, and a denoting a residual gauge freedom within harmonic gauge. 
Matching inner and outer metrics at r = a in the standard way yields 

m = M[l + (M/a) 2 + 4 (M/a) 4 + 0((M/a) 5 )}, (3a) 
a = 6[1 + (M/a) + (M/a) 2 + 2 (M/a) 3 + 0((M/ a) 4 )], (3b) 

P = -^{l + (M/a)~2(M/a) 2 + 0((M/af)}, (3c) 

C - --(M/af--(M/af--L(M/af + 0((M/af)]. (3d) 
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This initial data contains a discontinuity in the second derivative of the metric at 
the matching surface r — a. After fully relativistic evolution, we find this matching 
discontinuity to be smoothed out if we use the maximal slicing condition or "K- 
driver" slicing condition 0, both of which are given by solving an elliptic equation 
for the lapse function. "K-driver" slicing smoothed out the matching discontinuities 
faster than the maximal slicing case. Other slicing conditions such as static lapse 
(harmonic gauge) and the algebraic slicing condition a = 1 + log 7 are shown to 
lead to high frequency noise caused by the discontinuity at the matching surface. 

This study suggests that the matching surface is not necessarily fatal, provided 
suitable slicing conditions are used. Such matching surfaces are expected to exist 
in constructing initial data using the PN formulation for full GR evolution. Results 
from simulations using such a PN initial data will be reported in the near future. 

We thank Mike Pati and Nils Andersson for discussions. This work was sup- 
ported in part by NSF PHY96-00507 and PHY 96-00049, and by NASA ESS/HPCC 
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